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Abstract 



O 

We present a critical reappraisal of the available results on the broken phase 



of A(<& 4 )4 theory, as obtained from rigorous formal analyses and from lat- 
tice calculations. All the existing evidence is compatible with Spontaneous 
Symmetry Breaking but dictates a trivially free shifted field that becomes 

a.: 

controlled by a quadratic hamiltonian in the continuum limit. As recently 



pointed out, this implies that the simple one-loop effective potential should 
become effectively exact. Moreover, the usual naive assumption that the Higgs 
mass-squared m\ is proportional to its "renormalized" self-coupling \r is not 
valid outside perturbation theory: the appropriate continuum limit has 
finite and vanishing A#. A Monte Carlo lattice computation of the A($ 4 )4 
effective potential, both in the single-component and in the 0(2)-symmetric 
cases, is shown to agree very well with the one-loop prediction. Moreover, its 
perturbative leading-log improvement (based on the concept of Xr) fails to 
reproduce the Monte Carlo data. These results, while supporting in a new 
fashion the peculiar "triviality" of the A($ 4 )4 theory, also imply that, outside 
perturbation theory, the magnitude of the Higgs mass does not give a measure 
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of the observable interactions in the scalar sector of the standard model. 
1. Introduction 

The so called effective potential of a quantum field theory is a powerful tool to get 
informations about the structure of its vacuum, whose analysis, in turn, gives insight into 
its physical content. A reliable determination of the effective potential requires a critical 
control of details not directly accessible to a perturbative calculation. A convenient approach 
to uncovering the non perturbative features of such details may start by trying to determine 
the continuum limit of some discrete version of the theory on a lattice. In this approach 
one exploits the correspondence of an Euclidean (imaginary-time) quantum field theory to 
a relativistic (real-time) theory, the Schwinger's functions of the former determining the 
Green's functions of the latter. While substantial experience has been already gained in the 
numerical analysis of lattice field theories, the reliability of statements about their continuum 
limits is still open to questions, due to the fact that the equivalence of different procedures to 
get those limits is not yet controlled by theoretically well defined conditions. This amounts 
to a shift of the problem from the domain of bare computation, including the evaluation of 
numerical errors and/or approximations, to that of interpretation, exploiting the connections 
of the numerical results with the formal theory or with suitable models of the continuum 
limit. 

The continuum A($ 4 )d theories of a self-interacting scalar field $ in D spacetime di- 
mensions are perhaps the simplest examples to study as limits of their versions on a lattice. 
Their physical interest stems from their role in the Standard Model, as Spontaneous Sym- 
metry Breaking (SSB) and the Higgs phenomenon |]J determine the masses of the vector 
bosons and some other essential features of the Electroweak Sector of the model 0, e.g. 
with <3> being a (non hermitian) isospin doublet. 

The scalar A($ 4 )£> theories are characterized by a rather peculiar D-dependence of their 
dynamical content: they have been proved to be "trivial" ( i.e. equivalent to a generalized 
free field theory) for D > 4, but not so ( i.e. truly interacting) for D < 4 [0-0. For D = 4 the 
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situation is a bit less definite, since, instead of a proof, one has several "evidences" (based, 
of course, on different methods of analysis, whose assumed correctness points out a lack of 
equivalence) suggesting, in some sense, either "triviality" or "nontriviality" , the latter being, 
admittedly, an almost ignored option 0. 

The commonly accepted evidence about the "triviality" of A($ 4 ) 4 is frequently assumed 
to support the view that the scalar sector of the standard model, dynamically generating 
the vector boson masses through the Higgs mechanism Q, can only be an effective theory, 
valid up to some cutoff scale. Without a cutoff, the argument goes, there would be no scalar 
self- interactions, and without them there would be no symmetry breaking f|||. This point 
of view also leads to upper bounds on the Higgs mass 



Recently [|Tl|,|T^], it has been pointed out, on the basis of very general arguments, that 
SSB is not incompatible with "triviality", at least in the case of the scalar field theory with 
quartic self-interactions. One could have a non-zero vacuum expectation value (VEV) of 
the field, (<&), yet find only non-interacting, quasi-particle excitations above the vacuum. 
This picture, while providing a consistent model-theory characterized by the coexistence of 
"triviality" and of a non-trivial effective potential P~3"|-pT|], also clarifies the meaning of some 
explicit studies of "triviality" in the broken phase ((<&) 7^ 0) P^| . 

As discussed in |TT|Jl^| , since the effective potential determines the zero-momentum 1PI 



vertices of the symmetric phase p5| , its structure might carry some traces of non trivial 
interactions modifying the zero-momentum-mode sector of the symmetric phase in a fashion 
which somehow de-stabilizes the symmetric vacuum and gives rise to SSB. 



This remark |TT,12|, immediately brings out the relevance of studying the massless A($ 4 )4 



theories, as, for them, the zero-momentum (p^ = 0) is a physical on-shell point (inaccessible 
to perturbation theory). 

If, in the continuum limit, all interaction effects of the symmetric phase can be reduced 
into a set of Green's functions expressible in terms of just the first two moments of a Gaussian 
(functional) distribution, then one gets a description of SSB consistent with "triviality" 
which does not contradict any analytical or numerical result Hl3|-^l|,f7fl . In this case the 



shifted field h(x) = — ($) becomes governed by a quadratic Hamiltonian and the 

simple "one-loop" potential turns out to be effectively exact. More precisely, the one-loop 
potential and the gaussian effective potential determine the same relation 

ml = 8tc 2 v 2 r (f ) 

between the physical Higgs mass and the physical field VEV vr. This turns out to 
be true as well in the post-gaussian approximations (see, in particular, f27j), where the 
Higgs propagator G(x,y) is properly optimized at each value of = (<&), by solving the 



corresponding non-perturbative gap-equation G = G Q ((f)) [ p!f| , p^ ,^7 Hi^ , |2l| , p7| ■ This result 
follows from the exact connection of the effective potential Ks (0) with the effective potential 
for composite operators E[4>, G] p8f 



Jd 3 xV cS (<f ) ) = E[ ( f )l G 



and from "triviality" , which entails that all interaction effects (apart from those included in 
the physical vacuum structure) should be reabsorbed into the physical mass of the "gener- 
alized free field" h(x) defined above. In any approximation consistent with the (assumed) 
exact "triviality" of the theory the mass of the Higgs field cannot represent a measure of the 
observable interactions in the broken phase which, in the absence of interactions involving 
other fields, should contain only trivially free excitations. 

As we shall explain in sect. 2, this simple picture of an interacting zero- momentum 
mode that produces a non-trivial shape for the effective potential, giving SSB, co-existing 
with trivially non-interacting excitations in the broken phase, is not in contradiction with 
any knowledge derived from both all formally rigorous theorems about A($ 4 )4 in the broken 
phase and all available lattice calculations. 

It is well kown, since a long time, that the effective potential of A($ 4 )£> can be analysed 
in a rather simple way on a lattice by computing the VEV of the field as a function of a 
(global or uniform) "external source" strength J P^I23|J2^ . 



We report, in sects. 3,4, an analysis of the data obtained from precise Monte Carlo lattice 



calculations of the effective potential for both the single component and the 0(2)-symmetric 
A($ 4 ) 4 theory. 

We have performed such an analysis for several values of the bare quadratic (r ) and 
quartic (A ) couplings in the Euclidean discrete action in four dimensions 

where x stands for a general lattice site and, unless otherwise stated, lattice units are 
understood. 

The numerical results from a version of Metropolis' Monte Carlo algorithm have been 
best fitted with analytical forms possibly consistent with "triviality" and allowing for some 
suitably chosen "corrections". 

As a matter of fact, the one- loop potential agrees remarkably well with the lattice results, 
while its perturbative "improvement" fails to reproduce the Monte Carlo data. 

Thus, definite support is obtained for the picture proposed in 0,0] . The basic physical 
interest of our analysis, emphasizing a rather unconventional but theoretically not surpris- 
ing peculiarity of "renormalization" in this Bose field example, also stems from the new 
interpretation given to the Higgs mass, no more appearing as proportional to a renormal- 
ized coupling constant (whose vanishing is required by "triviality"), with remarkable new 
implications on the dynamics of the Higgs particle in the Standard Model. 

As we shall explain in sect. 2, our numerical results agree with those already obtained 
from lattice calculations, as far as we know, and our theoretical analysis is fully consistent 
with the general conditions imposed on the theory by the coexistence of SSB and triviality. 

Our analysis also shows that "triviality" and perturbation theory, within any attempt to 
construct a cutoff-independent and non-vanishing renormalized coupling at non zero external 
momenta Xr, are in inner contradiction ||29|| . The (new) evidence of this incompatibility is 
given by the lattice computation of the effective potential which turns out to represent 
just the sum of the classical potential and of the zero-point energy of a free field. Outside 
perturbation theory, and Xr are not proportional and while the latter vanishes the 
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former remains an arbitrary parameter. This result has important physical implications as 
discussed in sect. 5 where we shall also propose our conclusions. 



2. SSB and "triviality" 

In our study of SSB the basic quantity to compute is the VEV of the bare scalar field 

= MJ) (3) 

in the presence of an "external source" whose strength J(x) = J in Eq.(2) is ^-independent. 

In fact, determining 4>b{J) at several J- values is equivalent p3[^4| to inverting the rela- 
tion 

involving the effective potential I4//(</>b). Actually, starting from the action in Eq.(2), the 
effective potential of the theory can be rigorously defined, up to an arbitrary integration 
constant, through the above procedure. 

In this framework, the occurrence of SSB is determined by exploring (for J ^ 0) the 
properties of the function 

MJ) = -M-J) (5) 

in connection with its behaviour in the limit of vanishing J 

lim MJ) = ±vb^0 (6) 

over a suitable range of the bare parameters r a , X appearing in the lattice action Eq.(2). 

The existence of the A($ 4 )4 critical point can be stated, in the framework connecting 
the lattice theory with its continuum limit, on the basis of rigorous formal arguments [|J. 
Namely, for any A Q > 0, a critical value r c = r c (A D ) exists, such that, for r Q < r c , one finds 
non trivial solutions of Eq.(6). After fixing the vacuum and hence ($) = vb, one can define 
a shifted field 



h(x) = — vb 



(7) 



whose VEV vanishes by construction. 

As reviewed in the last section of ||, the following statements hold true. For D > 4, 
the continuum limit, precisely defined in Ch.15 of ||, is "trivial" in the sense that the 
connected n-point Green's functions (n > 3) vanish in that limit. Moreover, the shifted field 
renormalization constant Zh, (as derived from the renormalization of the term containing the 
field derivatives in the bare action) must be finite. As a consequence, if one assumes a single 
renormalization constant for both the field VEV vb and the shifted field h(x), one cannot 
get both a finite Higgs mass and a finite vacuum expectation value in the continuum limit 
||. Indeed, the bare field VEV vb becomes infinite in units of the mass of the shifted 
field which, in this specific context, must be trivially free. 

In this picture, where the same U Z" controls the scaling properties of both vacuum field 
and fluctuations, namely 



(8) 



what amounts to agree that 



the above quoted statements pfl requiring 



mi 



Z h 



(9) 



(10) 



and the finiteness of Zh = Z, also entail 



m l Z h 



mt 



(11) 
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In this approach, where ~ v \ is kept fixed and related to the Fermi constant, m\ vanishes 
in the continuum limit and one concludes that the theory makes sense only in the presence 
of an ultraviolet cutoff. This point of view, also leads to upper bounds on the Higgs mass 
|. In all interpretations of the lattice simulations performed so far, the validity of 



Eq.(8) has been assumed. Let us refer to the very complete review by C.B.Lang ||. There, 
the value of Z = is extracted from the large distance decay of the Goldstone boson 
propagator and used in Eq.(9) to define a renormalized field VEV vr from the bare field 
VEV vb as measured on the lattice. Since the lattice data both support the trend in Eq.(lO) 
and provide Z h = 1 to very good accuracy (see Tab. II of 0) one deduces (11). The limiting 
behaviour (11), obtained in this way, is usually interpreted (see |§ Fig. 19, caption included) 
to be a test of the validity of the perturbative relation 

m 2 h = \ R v% (12) 

in which the renormalized coupling A#, at external momenta comparable to the Higgs mass 
itself, has been introduced. 

As discussed in the Introduction, X R is a spurious concept in a "trivial" theory. On the 
other hand, without assuming the validity of any perturbative or quasi-perturbative relation, 
the lattice data confirm the conclusions of ||, namely, in the continuum limit 

a) the bare field becomes infinite in units of the Higgs mass, i.e. (10) 

b) the shifted field (fields) is (are) trivially free since Z^ — > 1 

Point a), supported by theoretical arguments and by the results of lattice simulations, 
can be taken as a basic condition defining the continuum limit of SSB in $4 theories. 

Point b) is completely consistent with the usual treatment, based on the Lehmann spec- 
tral representation of the propagator, for a field with vanishing VEV. To this end one 



identifies |3(J the Zh of the continuum theory, related to the (integral of the) spectral func- 
tion p(s), with the limit- value of its cutoff-dependent form when the ultraviolet regulator is 
(continuously) removed. Thus Zh = 1 implies p(s) to be just a ^-function, since the relevant 
fields have a vanishing VEV. 

Points a) and b) describe two basic outcomes of "triviality" in the broken phase. A 
consistent description of SSB in the A(<£> 4 )4 theory should explain the non-uniform scaling 
of the Higgs mass and of the bare vacuum field while preserving, at the same time, the non 
interacting nature of the shifted field(s). 
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Does perturbation theory provide any clue to understand a) and b) ? Strictly speaking, 
the perturbative one-loop /^-function, exhibiting an unphysical Landau pole, does not define 
any continuum limit. The usually accepted view maintains that the interval between the 
origin and the Landau pole becomes arbitrarily large in the limit A.r — ► 0, thus allowing to 
recover the continuum limit of the theory in agreement with "triviality" while still preserv- 
ing, at finite ultraviolet cutoff, the validity of the two perturbative relations (9,12). This 
interpretation of "triviality" can be questioned for the following reasons. In perturbation 
theory, Xr is the basic Renormalization Group (RG) invariant whose cutoff-independence 
generates the cutoff-dependence of the bare coupling A D = A Q (a), with a denoting the lat- 
tice spacing equivalent to the ultraviolet cutoff A. Different A^'s define different theories. 
It would be contradictory, and inconsistent with the original premise that Xr is finite and 
cutoff-independent, to now conclude that Xr has any implicit or explicit cutoff dependence 
that causes it to vanish in the limit a — > 0. "Triviality" means that Xr, which represents a 
measure of the observable interactions, must vanish identically in order to get a flow of the 
bare parameters r a = r Q (a), X Q = A c (a) extending down to a = 0. The one-loop perturbative 
/3-function, obtained under the assumption Xr > (and ^% In ^ < 1 ) simply does not 

h 

allow for such an extension. Furthermore, the above perturbative interpretation of "trivi- 
ality" is inconsistent with explicit two-loop calculations, where the perturbative /5-function 
exhibits a non trivial ultraviolet fixed point at non zero bare coupling A Q = A* > (let 
us remark that rigorous arguments || can be given to qualify as "spurious" this additional 
fixed point). In the perturbative approach, if Xr is taken to represent a value of the bare 
coupling at some finite momentum scale, A G = A D (a) flows towards A* when a — > for any 
< \ R < A*. 

Moreover, the perturbative explanation of Eq.(10), namely m\ — > at v B ~ |p = fixed, 
is certainly not unique. One might as well consider, e.g. the alternative possibility v 2 B — > 
+oo, m\ = fixed, if the renormalized vacuum field vr and vb are non trivially related 
through a renormalization constant ^ Z^. This alternative, first proposed in |TJ and 



later on thoroughly discussed in |11,12|J2y], does not violate any basic principle. Let us 



emphasize that the vacuum field renormalization Z§ has nothing to do, in principle, with 
the cutoff-dependence of the term containing the field derivatives in the bare action. The 
usual assumption of a single "field" -renormalization factor derives from giving an operatorial 
meaning to the field rescaling, namely 

This relation is a consistent shorthand for expressing the "wave function" -renormalization 
in a theory allowing an asymptotic Fock representation, but no more so in the presence of 
SSB, since, then, it overlooks that the shifted field is not defined before fixing the vacuum. 
The Lehmann spectral-decomposition argument constrains only Z^. 

The previous discussion should make clear that the perturbative interpretation (9,12) of 
the basic results a) and b) might be hardly taken as correct. At the same time, it should 
make apparent that the problems associated with the peculiar structure of the continuum 
limit of the A($ 4 )4 theory must be faced anew with open minded critical attention. 

In the following we shall briefly review the nonperturbative renormalization pattern 
which arises from the analysis of the effective potential in those approximations which are 
consistent with "triviality", as referred to in the Introduction and explained in |TT|,|T2" 



This will clearly show how the two basic results a) and b) are consistently derived in our 
approach. The basic difference with respect to perturbation theory, where one attempts to 
define a continuum limit in the presence of observable interactions, lies in the choice of the 
quantities that are required to remain finite in the continuum limit: namely, the vacuum 
energy and the "gap" in the energy spectrum corresponding to the mass of the /ree-field 
quanta. 

In terms of the bare vacuum field <Pb = and of the bare coupling A G we obtain the 
well-known result for the one-loop effective potential |25[ 



^-'~^H>U^(m^-i) (13) 

where, by definition, uj 2 (4>b) — 3A O 0^ and A denotes the Euclidean ultraviolet cutoff equiv- 
alent to the lattice spacing (A ~ -). The minimum condition requiring 

10 



m\ = uj 2 {vb) = 3X Vg = A 2 exp( 



16tt 2 
"9A7 



(14) 



with = |^-, also allows giving the one- loop effective potential in the form 



V 1 " 1oop (0b) 
Finally, the ground state energy is 



7T 



2^4 



'B /, VB - n 
( ln 72-- o)> 



1 



Z 2 v v?, 2 



(15) 



1287T 2 



(16) 



By defining the physical vacuum field 4>r — ~z~ i n such a way that all zero-momentum 



interaction effects are reabsorbed into its normalization, i.e. so that 

^yl-loop 



^2 \4>R = ±VR 



(17) 



holds true, we get P,P,0,|mp 



V 1 - toop (^) = 7r 2 ^(ln^-i 



(18) 



and 



mi 



-22 



57T U 



(19) 



The non-perturbative nature of the vacuum field renormalization {Z$ ~ 1/ A D ) , first discov- 
ered in the gaussian approximation by Stevenson and Tarrach should not be confused 
with the /i-field wave function renormalization. In fact, at one loop, h is just a free field 
with mass uj(<J)b), arid hence one has trivially Z^ = 1, in agreement with the basic result b). 
The structure with Z^ ^ Zh, is allowed by the Lorentz-invariant nature of the field decom- 
position into p^ = and p^ ^ components HTT] , |T~2|| . It is, of course, more general than the 
one assumed in perturbation theory, and its unconventional appearance is, in a sense, an 
indication of the peculiar non perturbative nature of the physical vacuum in the massless 
A($ 4 )4 theory, which allows SSB to coexist with "triviality". Notice that the possibility of 
a different rescaling for the vacuum field and the fluctuations is somewhat implicit in the 
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conclusions of the authors of ( see their footnote at page 401: "This is reminiscent of the 
standard procedure in the central limit theorem for independent random variables with a 
nonzero mean: we must subtract a mean of order n before applying the rescaling n~ l l 2 to 
the fluctuation fields"). 

To further illustrate the asymmetric treatment of the zero-momentum mode ( which is 
well known to play a special role in a Bose system at zero temperature) one can consider the 



analogy with the quantum mechanical 3-dimensional 5- potential discussed by Huang f24 



and Jackiw [ 3l| . In this example, the exact solution of the Schrodinger equation requires 
the radial wave-function to vanish at the origin. This condition is automatically satisfied 
for all partial waves except S-waves. For S-waves it cannot be satisfied if the wave function 
is required to be continuous at the origin. In that case, there would be no S-waves, and the 
solutions of the equation would not form a complete set. But a discontinuity at the origin 
might be acceptable, since the potential is singular there. By defining the ^-potential as the 
limit of a sequence of square well potentials, it is straightforward to show that, in the limit 
of zero-range, the S-wave solution vanishes discontinuously at the origin. Not surprisingly, 
all scattering phase shifts vanish. The analogy with A$ 4 theory is remarkable. There is 
no observable scattering, and yet the particle is not completely free. It is free everywhere 
except at the origin where the S-wave vanishes abruptly. In our field theoretical case, all 
finite momentum modes are free and as such there is no non-trivial S-matrix. However, 
the unobservable zero-momentum mode is a truly dynamical entity whose non-trivial self- 
interaction is not in contradiction with "triviality" being reflected in the zero-point energy 
of the free shifted field. 

In conclusion, the basic result a) is in complete agreement with the analytical results de- 
rived from the one-loop effective potential. Indeed, the finiteness of the RG-invariant ground 
state energy in Eq.(16) requires that the Higgs mass has also to be finite. In the continuum 
limit one gets a cutoff-independent rrih and a finite ground state energy in connection with 
the vanishing of the bare coupling as 
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In this limit, from Eq.(14), one rediscovers, with Z h = 1 identically, the trend in (10) in 
complete agreement with the previous points a) and b). Eqs.(18.19) also hold in the more 
sophisticated gaussian and postgaussian approximations [|TT|,0,0,|I^,^j . It is important to 
note that the logarithmic trend in Eq.(20), controlling the cutoff dependence of the ratio 

V B 

(see Eq.(14) ), allows to explain the observed logarithmic decrease on the lattice of the ratio 

m\Z h m\ 



v% v\ In K 

m h 

[see fig. 19 of f9j) without any need of introducing the leading-log Xr 

A A " 



1 + 9^ ln 

h 



Indeed, the usual perturbative interpretation of "triviality" (see fig.6 of || ), based on the 
assumed reliability of the above equation at A Q = +oo, i.e. 

16tt 2 



A* 



91n4' 



neglects that the leading-log resummation is not defined in that region and, as previously 
recalled, contradicts explicit perturbative two-loop calculations. 

3. Lattice computation of 0s(J). 

The VEV of the bare field = 4>b(J) is the simplest quantity to compute on the 
lattice, and its J-dependence can be exploited to get the first derivative of the effective 
potential, Eq.(4), as already recalled at the beginning of sect. 2. In this sense, a computation 
of the slope of the effective potential, based on the response of the system to an applied 
external source, has the advantage of being fully model-independent. At the same time, 
since it is well known P|,|3"2"|| that the data become unstable for small values of the external 



source, we have to restrict our analysis to a "safe" region of J (in our case \J\ > 0.05) where 
finite size effects appear to be negligible as checked on several lattices (see below). Indeed, 
at smaller |J| the precision of the data becomes uncertain (not only but especially in the 
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neighborhood of the critical point) and this was confirmed by testing the accuracy of the 
data against Eq.(5). The loss of accuracy observed in the neighborhood of J = cannot be 
accounted for by numerical and statistical errors but is sensitive to finite size effects. In this 
way we have been led to avoid any "direct" lattice calculation of <Pb{J) at J = 0, at least 
in a "small sized" lattice as 10 4 , and this is the main motivation for our choice of a method 
not relying on such "direct" calculations in the analysis of the lattice data. 

In order to compare the lattice data with the continuum theory one has to employ some 
definite model. As recalled in the Introduction, the equivalence of different procedures to get 
this limit is not yet safely controlled by theoretically well defined conditions. For instance, 
as discussed by Lang one can use an effective theory (the so called "chiral perturbation 
theory" ) to relate the lattice quantities computed in a finite volume at non-zero J to those 
at J = in an infinite volume. The cutoff-dependence of such "thermodynamic" quantities 
can be then compared with models of the continuum limit. 

In our case, we shall use two basically different functional forms, based on the one-loop 
effective potential and its perturbative leading- log improvement, to fit the measured values of 
4>b{J) and to extrapolate <Pb{J) toward J — > 0^. Our choice is motivated by the discussion 
of SSB and "triviality" presented in sect. 2. However, anybody could use our data, collected 
in a "safe" region of J and determined in a model- independent way, to compare with his 
own preferred model and draw the appropriate conclusions. 

The calculation of 0#( J), on the basis of the action given in Eq.(2), has been performed by 
using a version of the Metropolis' Monte Carlo algorithm with periodic boundary conditions, 
for several values of the parameters r a , A Q , J appearing in the discrete action. Other currently 
used algorithms, as the cluster method and the hybrid Langevin's approach, are presently 
under our investigation and the results will be compared in a forthcoming paper. However, 
possible improvements are only expected in the region of very small J which is not analyzed 
here. 

The calculations have been carried out on DEC 3000 AXP (Mod. 400 and 500) computers. 
With the aim of obtaining a reliable control on the precision of the data, we have carefully 
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analyzed their dependence on the choice of the random number generator, of the "thermal- 
ization" process and of the dimension of the sample (number of "sweeps" or "iterations") 
entering the determination of the relevant data. The results for J = J(4>b) obtained from 
runs on our 10 4 lattice have been compared with those reported in p^j2^| and the agreement 
with the plots J = J(4>b) at A = 1 in Fig.2 of [F23J (suitably enlarged) has been found better 
than 1% in the range \J\ > 0.05. Other comparisons with [p3jp4| have not been attempted 



in view of their direct determination of vb at J = and of their assumption = Zh in the 
evaluation of the shifted field propagator. 

After 60,000 iterations, which correspond to the data reported in Table I, the calculated 
4>b{J) is stable at the level of the first three significant digits for all | J |- values. This stability 
has been controlled with several runs up to 200,000 iterations. As | J\ increases, the stability 
gets better (down to the level of 10~ 4 in the explored range of J-values) and finite size effects 
appear negligible, as checked by comparing with the runs on a 16 4 lattice and, subsequently, 
confirmed by an independent computation performed by P. Cea and L. Cosmai |$3| reported 
below. 

In order to introduce our analysis of the Monte Carlo data, let us first discuss the 
predictions from the one- loop potential which, as discussed in sect. 2, reflect the "triviality" 
of the theory. By differentiating Eq.(15), we obtain the "bare source" 



J{*b) = - -T7 ^ =^ln^, (21) 



which can be directly compared with the lattice results for J = J(0#). Strictly speaking, 
just as the effective potential is the convex envelope of Eq.(15) |34l , |2T|l , Eq.(21) is valid only 
for \</>b\ > vb and J(4>b) should vanish in the presence of SSB in the range — vb < 4>b < vb- 
This means that the average bare field 4>b{J) — (<&b)j, while satisfying the relation (5), i.e. 

MJ) = -M-J) ( 5 ) 

for any J ^ 0, should tend to the limits stated in Eq.(6). However, as discussed above, 
Eq.(5) is not well reproduced numerically for \J\ ~ 0.01 or smaller and we shall have to 
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restrict our analysis to a "safe" region of J- values ( | J\ > 0.05) allowing to get the values of 
v B and from a best fit to the lattice data by using Eq.(21) after having identified on the 
lattice the "massless" regime, i.e. the one corresponding to a renormalized theory with no 
intrinsic scale in its symmetric phase ($) = 0. 

Usually this would require to find the value of the quadratic bare coupling r a for which 

n 



7T\4>b=0 - 



d 2 V 

However, the region around <pB — being not directly accessible, we have argued as follows. 



We start with the general expression |T2 



J(cf> B ) = a<f>% H<j> 2 B /v B ) + PvlfaO- ~ 01/4), (22) 

which is still consistent with "triviality" (corresponding to an effective potential given by 
the sum of a classical background and the zero point energy of a free field) but allows for an 
explicit scale-breaking /9-term. Setting a = one obtains a good description of the data in 
the "extreme double well" limit (r D much more negative than r c , where r c corresponds to the 
onset of SSB) where SSB is a semi-classical phenomenon and the zero-point energy represents 
a small perturbation. Then we start to increase r Q , at fixed A Q , toward the unknown value r c 
and look at the quality of the fit with (o£,/3,vb) as free parameters. The minimum allowed 
value of r such that the quality of the 2-parameter fit (a, (3 = 0,v B ) is exactly the same 
as in the more general 3-parameter case will define the "massless" case so that we can fit 
the data for J = J{4>b) with our Eq.(21). Finally, Eq.(14) suggests that the vacuum field 
Vb, in lattice units, vanishes extremely fast when A Q — > 0. Thus, to avoid that noise and 
signal become comparable, we cannot run the lattice simulation at very small values of A G 
but have to restrict to values A G ~ 1 such that still % << 1 but v B , in lattice units, is 
not too small. Smaller values of A G (~0.4-0.6), however, should become accessible with the 
largest lattices available today (~ 100 4 ) where one should safely reach values \J\ ~0.001 
being still in agreement with Eq.(5). At A D = 1 we have identified the massless regime at a 
value r = r s < r c where r s ~ —0.45. By using the accurate weak-coupling relation between 
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the bare mass and the euclidean cutoff fl25 



3Ao r A 2 (23) 



167T 2 

and using the relation A = ixyq (where yq is expected to be 0(1)) we obtain yq ~ 1.55. It 
should be noted that our operative definition of the "massless" regime amounts to determine 
for weak bare coupling the numerical coefficient relating ultraviolet cutoff and lattice spacing 
for quadratic divergences (A ~ 4.87 ). This agrees well with an independent analysis of the 
Liischer and Weisz lattice data presented by Brahm ||35|| . In fact Brahm's result (in the 
range A < 10 ) is A = 4.893 ± 0.003. Also, ref. [3l| predicts the massless regime to 
correspond to r s = —(0.224 ± 0.001) for A Q = 0.5 in the infinite-volume limit. This implies 
r s ~ —(0.448 ± 0.002) for A G = 1, in excellent agreement with our result r s ~ —0.45. Thus, 
the model-dependence in the definition of the massless regime, introduced by our fitting 
procedure to Eq.(22), is negligible for ^2 << 1- 
Finally, the massless relation (14) predicts, 

K) 1 - loop = ^exp(-^) (24) 



/3X 9A 
and 

Z 1 ^ = |£. (25) 

In Eq.(24) we have used i/l rather than yq since one does not expect precisely the same coef- 
ficient to govern the relation between euclidean cutoff and lattice spacing for both quadratic 
and logarithmic divergences. 

The consistency of Eqs. (24,25) with the corresponding values obtained from the best 
fit to the lattice data checks the validity of the picture given in sect. 2 for the coexistence 
of SSB and "triviality". We have determined y^ at A = 1 and, after that, the scaling 
behaviour of vb with the bare coupling has been directly compared with Eq.(24). On the 
other hand, the prediction in Eq.(25) for Zq, which fixes the absolute normalization of the 
curve J = J(0s), only depends on A D and can be directly compared with the lattice data in 
the massless regime. 

17 



Our numerical values for <Pb{J) and the results of the 2-parameter fits to the data by 
using Eq.(21) for A Q = 0.8,1.0 and 1.2 are shown in Table I together with the one-loop 
prediction (25). The change of r Q with A D is computed by using Eq.(23) and our result 
r s a 2 ~ —0.45 for A Q = 1. The reported errors are due to three main sources : 

a) the purely statistical errors from the Monte Carlo sampling of the field configurations 

b) errors which can be easily recognized as numerical artifacts (in our case they can 
be extimated by comparing with Eq.(5) ) 

c) finite size effects which we have extimated by computing a few points on a 16 4 

lattice 

By fixing to its one-loop value (25) and using Eq.(24) for vb we have determined the 
value y L = 2.07±0.01 ( -^j = which corresponds to v B = (5.814±0.028)10~ 4 at A = 1. 
Once we know we can predict the value of vb = {vb) T)i and compare with the results of 
the fit to the lattice data at A o =0.8 and 1.2 . The results in Table II show a remarkable 
agreement of the Monte Carlo data with the one-loop predictions. 

As an additional check that our values for 0_b(J) are not affected by appreciable sys- 
tematic effects, we have compared with a completely independent Monte Carlo simulation 
performed by P. Cea and L. Cosmai [33|| on a 16 4 lattice. Their data for A D = 1 and 



r Q = —0.45 are reported in Table III. Note that their errors are statistical only and are 
smaller by an order of magnitude than those reported in our Table I (where errors of type 
b) and an estimate of the errors of type c) were included). Still, the quality of the one- loop 
fit is excellent confirming that finite size effects do not affects the results of the fits. 

At this point one may ask: " How does the standard perturbative approach compare with 
the lattice data ?". This is not a trivial question. Indeed, the same lattice data might be 
consistent with different functional forms thus allowing for different extrapolations at J = 
and, hence, different interpretations. In particular, since our data are collected in the region 
of small bare coupling (^f << 1 ), it is reasonable to check whether the lattice data agree 
better with the one-loop potential or with its perturbative leading-log improvement, related 
to the concept of \r. This is even more important in our case because: 



i) The leading-log approximation provides a commonly accepted perturbative inter- 
pretation of "triviality" (ignoring the objections raised in sect. 2) 

ii) Just the validity of the leading-log resummation is usually considered as an indication 
that the one-loop minimum cannot be trusted and there is no SSB in the massless regime 
|25fl. If this were true the correct estrapolation at J = would require vb = 

However, one should be aware that the higher order corrections to the one-loop potential 
represent genuine self-interaction effects of the shifted field which we know from "triviality" 
cannot be physically present in the continuum limit. Were we observing the need for per- 
turbative corrections of the one-loop formula this would contradict all the non perturbative 
evidence about h{x) being a generalized free field. 

By fitting our lattice data to the leading-logarithmically "improved" version of Eq.(21), 
namely 



K4>b 

1 + Ik In » 



(xll denoting an adjustable parameter for each A G ) we find, for 21 degrees of freedom, 
(x 2 )ll = 53,163,365 for A Q = 0.8,1.0,1.2 respectively (to compare with the values 
X 2 =5,6,13 of the corresponding one-loop 1-parameter fits). The same comparison with 
the Cea-Cosmai data reported in Table III gives the result 

x LL = 1.627 ±0.001 

x 2 moo 

cLol ~ 9 - 1 

This should be compared with the one-loop results = reported in Table III for the 
2-parameter fit and the value -^j = JMj- for the 1-parameter fit when is constrained 
to its one-loop value = 26.319. These results show that the agreement between lattice 
data and one- loop formulas are not a trivial test of perturbation theory but, rather, a non 
perturbative test of "triviality" . They also verify the consistency of our definition of the 
continuum limit presented in sect. 2. 
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In order to fully appreciate what the above results exactly mean we have reported in 
Tables IV and V the two theoretical sets J = J Th (4>B) (Th=one-loop, leading-log) which best 
fit the measured values of 4>b (with their errors A J corresponding to the error in 4>b)- Finally, 
in the last column, we report the quantity x 2 = (</ EXP — J Th ) 2 /(AJ) 2 where J = J EXP is 
the "experimental", input value of J at which the various 4>b have been computed by Cea 
and Cosmai. By inspection of Tables IV and V one finds: 

1) the leading-log formula is accurate to the level of a few percent. The fit routine 
tries to reproduce the lattice data with a suitable value of the free parameter xll within 
the logarithm, thus optimizing the relation between euclidean cutoff and lattice spacing and 
effectively reabsorbing also non-leading perturbative corrections. In this case the experi- 
mental data and the leading- log curve intersect each other somewhere in between J =0.5 
and J = 0.4. However, their slopes are clearly different as it can be seen from the fact that 
the theoretical prediction is lower than the experimental data for J > 0.5 and higher for 
J < 0.4. Since the Monte Carlo data are so precise the value of the x 2 is enormous. 

2) the one-loop prediction is accurate to a completely different level, namely O(10 -4 ). 
This result cannot be understood from the conventional perturbation theory viewpoint; 
apparently a cruder approximation reproduces the lattice data remarkably well, while its 
"pro forma" improvement fails badly ! However, if one understands our interpretation of 
"triviality" , then this is precisely what one expects. 



In conclusion, the lattice computations nicely confirm the conjecture [|n],|ll|] that, as a 
consequence of "triviality" , the one-loop potential becomes effectively exact in the continuum 
limit. Indeed, there is a well defined region in the space of the bare parameters (r D , A D ), con- 
trolled by Eq.(23), where the effective potential is described by its one-loop approximation 
to very high accuracy. The lattice analysis of the effective potential in this region, i.e. close 
to the critical line and in the broken phase, where the theory is known to become gaussian, 
should be pursued with larger lattices and, possibly, also with different methods. Indeed, 
by increasing the lattice size, our analysis can be further extended towards the physically 
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relevant point ^§ — > + , r c — > which corresponds, in the continuum limit of the theory, 



to "dimensional transmutation" [25| from the classically scale invariant case. Our results 



confirm the point of view []Il] , fL2 | that the massless version of the A($ 4 )4 theory, although 



"trivial", is not "entirely trivial": it provides, at the same time, SSB and a meaningful 
continuum limit A — > oo, A„ — ► + such that the mass of the free shifted field in Eq.(14) is 
cutoff independent. However, the mass rrih does not represent a measure of any observable 
interaction since the corresponding field h(x) is described by a quadratic hamiltonian. 

4. The 0(2)-symmetric case 

As discussed in [|T7| , |T8|, [TT|1 , and explicitly shown in f27|, one expects Eqs. (18,19), of the 



one-component theory, to be also valid for the radial field in the 0(N)-symmetric case. This 
observation originates from ref. |36| which obtained, for the radial field, the same effective 
potential as in the one-component theory. This is extremely intuitive. The Goldstone-boson 
fields contribute to the effective potential only through their zero-point energy, that is an 
additional constant, since, according to "triviality", they are free massless fields. Thus, in 
the 0(2)-symmetric case, one may take the diagram {V e ff, 4>b) for the one-component theory 
and "rotate" it around the V e ff symmetry axis. This generates a three- dimensional diagram 
(V e ff, 0i, 2 ) where V e ff, owing to the 0(2) symmetry, only depends on the bare radial field, 



Pb = vV? + # (26) 

in exactly the same way as V e ff depends on <pB in the one-component theory; namely 
(u\p B ) = 3X oP l) 

T/ l-loop/- \ _ A ° 4 , ^(Pb) ( UJ 2 (p B ) 1\ 

By using Eqs. (14,25) ^/ 1_loop can be re-expressed in the form 

V 1 - loo »(PB) = ^(ln P 4- 1 -), (28) 



By differentiating Eq.(28), we obtain the bare "radial source" 
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j(pb) = ^M = ^A (29) 

dp B Aa, v B 



which we shall compare with the lattice results for J = J(pb)- 

The lattice simulation of the 0(2)-invariant theory has been obtained, with periodic 
boundary conditions, from the action 

E E(^(^ + ^) - $^)) 2 + E($, 2 W) + y (E W) 2 - E j,^ W] (so) 

x / Ai =lj=l A i=l 4 i=l i=l 

where $i and $2 are coupled to two constant external sources J\ and J2. By using J\ = 
Jcos9 and J 2 = J sin 9 it is straightforward to show that, having defined 0i = ($1)^,73 and 
02 = {^2)j 1 ,j 2 i t ne bare radial field 



p B = + 0| (31) 

does only depend on J, that is 

Pb = Pb(J) (32) 

We started our Monte Carlo simulation on a 10 4 lattice by investigating first the (r Q , A ) 
correlation corresponding to the classically scale-invariant case. 

In the 0(2) case, we have used the results of the previous section and of ]35| to get the 



analogous of Eq.(23) for the two-component theory. Since we are considering ^§ « 1, this 
simply introduces a combinatorial factor of 4/3 so that Eq.(23) becomes 



4tt 2 



A (33) 



Hence, we expect the massless case to correspond to r s ~ —0.6 for A G = 1. This was 
confirmed to good accuracy by using the fitting procedure to the analogous of Eq.(22) 
defined with ps in the place of 0b- Thus, the identification of the massless regime on the 
lattice appears to obey the simple scaling laws (23,33) and is under theoretical control. 

Just like in the single-component case Eq. (5) was poorly reproduced numerically at small 
values of J, in the 0(2)-symmetric case we find that, at small J, the exact ^-independence of 
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Pb (see Eq.(32) ) is poorly reproduced and a "direct" data processing becomes unreliable. We 
therefore consider a "safe" region of J-values, J > 0.05, in which the spurious ^-dependence 
gives the data an uncertainty less than ±3%. Fitting the data to Eq.(29) we can infer the 
values of vb and and compare with Eqs. (24,25). 

Our numerical values for ps(J) in the massless case, obtained with our version of the 
Metropolis' Monte Carlo algorithm on the 10 4 lattice, are reported in Table VI for A G =1.0, 
1.5 and 2.0. For each A G the corresponding r Q is computed by using Eq.(33) and our numerical 
input r s = —0.6 for A G =1 . The errors in Table VI are essentially determined by the observed 
spurious variation of ps in the range < 6 < 2n. As discussed above, this is a numerical 
artifact and should be considered a systematic effect of the Monte Carlo lattice simulation. 
It is reproduced with three different random number generators, all consistent with the 
Kolmogorov-Smirnov test at the level O(10~ 4 ). At low J this systematic effect completely 
dominates the error; the statistical errors, already after 30,000 iterations, are 4-5 times 
smaller. 

Table VI also reports the vb and values obtained from the two-parameter fits to 
the data using Eq.(29). The resulting Z^ values agree well with the one-loop prediction of 
Eq.(25). As in the one-component theory, to perform a more stringent test of the one- loop 
potential, we have next constrained to its one-loop value in Eq.(25) and, under such 
condition, we have done a precise determination of vb from a one-parameter fit to Eq.(29). 
This allows a meaningful comparison with Eq.(24). In this case the value of yi derived from 
the one-parameter fit to the 0(2) lattice data at A G = 1.5 and r = —0.9 turns out to be 
y L = 2.44 ±0.03 with -^j = jj. In Table VII we show the results of the one-parameter fits 
to the data at A G =1.0 and 2.0 and the comparison with Eq.(24) for y L = 2.44 ± 0.03. It is 
apparent from Table VII that the one-loop potential reproduces quite well the scaling law 
of vb with the bare coupling, as previously discovered in the one-component case. 

The value of yi obtained from the lattice simulation of the 0(2)-symmetric massless A$ 4 
theory is ~17% larger than the value yL = 2.07 ± 0.01 obtained in the one-component case. 
This is partly due to our choice of fixing the value of r s = —0.6 at A G = 1 with the simple 

23 



combinatorial factor 4/3 discussed above and ignoring finite size corrections to Eq.(33), see 



With a two-component field, the errors (as estimated from the spurious ^-dependence 
of pb) are larger than in the one-component theory. The massless regime, identified on the 
basis of the (J,Pb) correlation, is found in a range of revalues around r a = —0.6 for A D = 1. 
Choosing instead r s ~ —0.585 would give yi ~ 2.07. 

We have analyzed the 0(2) lattice data in the same fashion as already described in the 
one-component case. Again, if the agreement with the one-loop formulas would be a trivial 
test of perturbation theory the data should agree at least as well, if not better, with the 
leading-log formula based on the perturbative /3-function 

J ll (pb) - Xop% 



\ _|_ } n TT^LL 

' It? p B 



(xll denoting an adjustable parameter to be optimized at each A ). However, when we fit 
the A G = 1.0, 1.5 and 2.0 data to this formula we find, respectively, for 17 degrees of freedom, 
(x 2 )ll =9, 44 and 133 (to compare with the values x 2 =0.8, 3.3 and 7.3 obtained from the 
1-loop one-parameter fits with Eq.(29) when is constrained to its value in Eq.(25)). Note 
that we are working in a region where the "A log" term is not small, being of order unity. 
Thus, the good agreement between the data and the one-loop formula cannot be due to the 
higher-order corrections (expected by perturbation theory) being negligibly small; it must 
be due to their complete absence. Without the "triviality" argument, this would be an 
incomprehensible miracle. 

Finally, to have an idea of the finite size effects, we have repeated the lattice simulation 
on a 16 4 lattice for A = 1 and r D = —0.6. The computation has been performed at the 
same values 9 =30, 60, 90, 120, 150, 180, 210, 240 and 270 degrees for the angular source 
and the data are reported in Table VIII. The spurious dependence of the radial field on 9 
is slightly smaller than in the 10 4 case and the agreement with the one loop prediction in 
Eq.(25) improves. 

In conclusion, the good agreement between lattice simulation and eqs. (24,25,29) provides 
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definite evidence that the dependence of the effective potential on the radial field in the 0(2)- 



symmetric case is remarkably consistent with the expectations pl|J12|| based on the analysis 
of the one-component theory and with the explicit analytical postgaussian calculation in 
p7| which numerically confirmed, to very high accuracy, the validity of Eq.(19) in the 0(2) 
and 0(4) case. 



5. Conclusions 

We have proposed a critical reappraisal of the available results on the broken phase of 
A($ 4 )4 theory, suggesting a consistent interpretation of the coexistence of "triviality" and 
SSB, whose validity is also confirmed by the lattice calculations. All the existing analytical 
and numerical evidences can be summarized in the two basic points a) and b) illustrated 
in sect. 2. They represent the essential ingredients for a consistent description of SSB in 
pure A($ 4 )4 theory. They are naturally discovered in those approximations to the effec- 
tive potential which are consistent with "triviality" (one-loop, gaussian effective potential, 
postgaussian approximations and so on) where the shifted field(s) appear as generalized free 
field(s). 

As discussed in sect. 2, the perturbative interpretation of the lattice data based on the 
leading-log approximation is by no means unique and suffers of internal theoretical draw- 
backs. Even more, in sects. 3 and 4 we have provided extensive numerical evidence that in 
the scaling region, clearly identified with our procedure, there is no trace of any Higgs-Higgs 
or Higgs-Goldstone interactions in the lattice data for J = J(4>b)- The measured effective 
potential, being remarkably consistent with the peculiar exponential decay in Eq.(24) and 
with Eq.(25), is just the sum of the classical potential and of the zero point energy of a free 
field. The corrections, expected on the basis of the leading-log approximation, are definitely 
ruled out by our results. 

The simple scaling laws associated with the one-loop potential allow to recover a "not- 
entirely trivial" continuum limit where, despite of the absence of observable interactions, 
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one still gets SSB as the basic ingredient to generate the vector boson masses in the standard 
model of electroweak interactions. The A($ 4 )4 theory is not "entirely trivial" and can be 
physically distinguished from a free field theory. Indeed it exhibits SSB and a first order 



phase transition [H],|5?H at a critical temperature T c depending on the vacuum energy in 
eq.(16). Also, the origin of SSB from a classically scale-invariant theory, allows to reconcile 
the results of ref. fTfl with "triviality" . 

Let us emphasize, now, that the mass of the Higgs particle does not represent, by itself, 
a measure of any interaction and the naive proportionality between and Xr is not valid 
outside perturbation theory. While the latter vanishes the former remains an arbitrary 
parameter. As discussed in |l^,ll],|12,|T8| substantial phenomenological implications follow. 



In our picture, which is consistent with all rigorous and numerical results available today, 
the Higgs particle can be very heavy and, still, be weakly interacting. In fact, in the 
absence of the gauge and Yukawa couplings, it would be trivially free. In the most appealing 
theoretical framework, where SSB is generated from a classically scale-invariant theory, we 
get m\ = 87r 2 w|j so that, by relating vr to the Fermi constant (modulo purely electroweak 
corrections which are small if the top quark mass is below 200 GeV [ |TT| , |rT| , |r8| ) we get 
m/j ~ 2 TeV. In general the Higgs can be lighter or heavier, depending on the size of the 
scale-breaking parameter in the classical potential |12|| . In any case, beyond perturbation 
theory, its mass is not related to its decay rate through X R . 

For instance, consider the Higgs decay width to W and Z bosons. The conventional 
calculation would give a huge width, of order Gpm| ~ for ~ 2 TeV. However, in 
a perturbatively renormalizable-gauge calculation of the imaginary part of the Higgs self- 
energy, this result comes from a diagram in which the Higgs supposedly couples strongly to a 
loop of Goldstone bosons with a physical strength proportional to its mass squared. On the 
other hand, if "triviality" is true, all interaction effects of the pure A$ 4 sector of the standard 
model have to be reabsorbed into two numbers, namely rrih and vr, and there are no residual 
interactions. Thus, beyond perturbation theory, that diagram is absent, leaving a width of 
order g 1 ni} l . Therefore, if "triviality" is true, the Higgs is a relatively narrow resonance, 
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decaying predominantly to tt quarks. Although the scalar sector must be treated non- 
perturbatively, one may continue to treat the gauge interactions using perturbation theory. 
Effectively, then, inclusive electroweak processes can be computed as usual, provided one 
uses a renormalizable gauge and sets the Higgs self-coupling and its coupling to the Higgs- 
Kibble ghosts (the would-be Goldstones) to zero. One should avoid the so-called "unitary 



gauge" where, as well known since a long time ||38|| , there is no smooth limit M w — > in 
perturbation theory, as a consequence of the longitudinal degrees of freedom. The possibility 
of decoupling the gauge sector from the scalar sector of the standard model in the limit 
(7 — > is crucially dependent on the fact that gauge invariance is carefully retained at any 
stage. By computing in a renormalizable gauge one never gets effects (such as cross- 
sections, decay rates,...) growing proportionally to the Higgs mass squared unless there is 
a coupling in the theory proportional to the Higgs mass at zero gauge coupling. Thus, the 
Higgs phenomenology, on the basis of gauge invariance, depends on the details of the pure 
scalar sector, namely "triviality" and the description of SSB. We believe that our picture, 
supported as we have shown by the lattice results for the effective potential, deserves a very 
serious attention. The conventional wisdom, built up on the perturbative assumption that 
the Higgs mass is a measure of its physical interactions, runs a serious risk of predicting a 
whole set of physical phenomena which exist only in perturbation theory. 

We thank K. Huang and P. M. Stevenson for many useful discussions and collaboration. 
We are indebted to P. Cea and L. Cosmai for communication of their preliminary numerical 
results. 



27 



REFERENCES 



[1] P. W. Higgs, Phys. Lett. 12, 132 (1964); Phys. Rev. Lett. 13, 508 (1964); F. Englert 
and R. Brout, ibid, 321; G. S. Guralnik, C. R. Hagen, and T. W. B. Kibble, ibid, 585; 
T. W. B. Kibble, Phys. Rev. 155, 1554 (1967). 

[2] S. Weinberg, Phys. Rev. Lett. 19, 1264 (1967) ; A. Salam, Proc. 8 th Nobel Symp., N. 
Svartholm ed. (Almqvist and Wicksell Stockholm, 1968) 367. 

[3] J. Frdhlich, Nucl. Phys. B200(FS4), 281 (1982); M. Aizenman, Phys. Rev. Lett. 47, 1 
(1981); A. Sokal, Ann. Inst. H. Poincare, 37, 317 (1982). 

[4] J. Glimm and A. Jaffe, Quantum Physics: A Functional Integral Point of View 
(Springer, New York, 1981, 2nd Ed. 1987). 

[5] R. Fernandez, J. Frdhlich, and A. D. Sokal, Random Walks, Critical Phenomena, and 
Triviality in Quantum Field Theory (Springer- Verlag, Berlin, 1992). 

[6] K. G. Wilson and J. Kogut, Phys. Rep. C12, 75 (1974); G. A. Baker and J. M. Kincaid, 
Phys. Rev. Lett. 42, 1431 (1979); B. Freedman, P. Smolensky and D. Weingarten, Phys. 
Lett. B113, 481 (1982); D. J. E. Callaway and R. Petronzio, Nucl. Phys. B240, 577 
(1984); I. A. Fox and I. G. Halliday, Phys. Lett. B 159, 148 (1985); C. B. Lang, Nucl. 
Phys. B 265, 630 (1986); M. Liischer and P. Weisz, Nucl. Phys. B 290, 25 (1987). 

[7] J. Pedersen, I. E. Segal and Z. Zhou, Nucl. Phys. B376, 129 (1992). 

[8] D. J. E. Callaway, Phys. Rep. 167, 241 (1988). 

[9] C. B. Lang, Computer Stochastics in Scalar Quantum Field Theory, to appear in 
Stochastic Analysis and Application in Physics, Proc. of the NATO ASI, Funchal, 
Madeira, Aug. 1993, ed. L. Streit, Kluwer Acad. Publishers, Dordrecht 1994. 

[10] J. Kuti, L. Lin, and Y. Shen, Phys. Rev. Lett. 61 (1988) 678; H. Neuberger, U. M. Heller, 
M. Klomfass, and P. Vranas, in Proceedings of the XXVIth International Conference 



28 



on High Energy Physics, Dallas, TX, August 1992. 

[11] M. Consoli and P. M. Stevenson, Resolution of the A$ 4 Puzzle and a 2 TeV Higgs Bo- 
son, Rice University preprint, DE-FG05-92ER40717-5, July 1993, submitted to Physical 
Review D. (hep-ph 9303256). 

[12] M. Consoli and P. M. Stevenson, Zeit. Phys. C63, 427 (1994). 

[13] M. Consoli and A. Ciancitto, Nucl. Phys. B254, 653 (1985). 

[14] P. M. Stevenson and R. Tarrach, Phys. Lett. B176, 436 (1986). 

[15] P. Castorina and M. Consoli, Phys. Lett. B235, 302 (1990); V. Branchina, P. Castorina, 
M. Consoli and D. Zappala, Phys. Rev. D42, 3587 (1990). 

[16] V. Branchina, P. Castorina, M. Consoli, and D. Zappala, Phys. Lett. B274, 404 (1992). 

[17] M. Consoli, in "Gauge Theories Past and Future - in Commemoration of the 60th 
birthday of M. Veltman", R. Akhoury, B. de Wit, P. van Nieuwenhuizen and H. Veltman 
Eds., World Scientific 1992, p. 81; M. Consoli, Phys. Lett. B 305, 78 (1993). 

[18] V. Branchina, M. Consoli and N. M. Stivala, Zeit. Phys. C57, 251 (1993). 

[19] R. Ibanez- Meier and P. M. Stevenson, Phys. Lett. B297, 144 (1992). 

[20] R. Ibanez-Meier, I. Stancu and P. M. Stevenson, "Gaussian Effective Potential for the 
U(l) Higgs model", Rice Preprint DOE/ER/05096-51, July 1992 (frep-ph 9207276|) . 

[21] U. Ritschel, Phys. Lett. B318, 617 (1993) 

[22] M. G. do Amaral and R. C. Shellard, Phys. Lett. B171, 285 (1986); M. Liischer and 
P. Weisz, Nucl. Phys. 295, 65 (1988); J. K. Kim and A. Patrascioiu, "Studying the 
continuum limit of the Ising model". University of Arizona preprint AZPH-TH/92-09, 
July 1992. 

[23] K. Huang, E. Manousakis, and J. Polonyi, Phys. Rev. D35, 3187 (1987); 



29 



[24] K. Huang, Int. J. Mod. Phys. A4, 1037 (1989); in Proceedings of the DPF Meeting, 
Storrs, CT, 1988. 



[25 
[26 

[27; 

[28 
[29 

[30 

[31 

[32 
[33 
[34 

[35 



[36 



S. Coleman and E. Weinberg, Phys. Rev. D7, 1888 (1973). 

L. Polley and U. Ritschel, Phys. Lett. B221, 44 (1989); R. Ibanez-Meier, A. Mattingly, 
U. Ritschel, and P. M. Stevenson, Phys. Rev. D45, 2893 (1992). 

U. Ritschel, Zeit. Phys. C63, 345 (1994) 

J. M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. D10, 2428 (1974) 

M. Consoli and P. M. Stevenson, "Triviality" and the perturbative expansion in A$ 4 
theory, Rice University preprint, March 1994, submitted to Physics Letters B. 

W. Zimmermann, in Lectures on Elementary Particles and Quantum Field Theory, 1970 
Brandeis Summer Institute, Vol.1, pag.395, M.I.T. Press, Cambridge MA, 1970. 

R. Jackiw, 5-Function Potentials In Two- And Three- Dimensional Quantum Mechanics, 
Contribution to the Beg's Memorial Volume, World Scientific Publ. Co., 1992. 

D.J.Callaway and D.J.Maloof, Phys. Rev. D27, 406 (1983). 

P. Cea and L. Cosmai, private communication. 

A. Dannenberg, Phys. Lett. B202, 110 (1980); V. Branchina, P. Castorina, and D. 
Zappala, Phys. Rev. D41, 1948 (1990). 



D. E. Brahm, The lattice cutoff for \<$>j and A$^, CMU-HEP94-10 preprint, [hep- 
lat/ 9403021| , March 1994. 



L. Dolan and R. Jackiw, Phys. Rev. D9, 3320, (1974). 



[37] G. A. Hajj and P. M. Stevenson, Phys. Rev. D 37, 413 (1988); N. M. Stivala, Nuovo 
Cimento A 106 (1993) 777. 

[38] A. Vainshtein and I. Khriplovich, Sov. Journ. Nucl. Phys. 13, 111 (1971). 



30 



J X = 0.8 r Q = -0.36 A G = 1.0 r Q = -0.45 A G = 1.2 r Q = -0.54 



-0.700 


-1.0024 ±0.0003 


-0.9401 ± 0.0003 


-0.8935 ±0.0003 


-0.600 


-0.9540 ± 0.0003 


-0.8950 ±0.0003 


-0.8512 ±0.0003 


-0.500 


-0.8997 ±0.0005 


-0.8444 ± 0.0003 


-0.8037 ±0.0003 


-0.400 


-0.8371 ±0.0005 


-0.7867 ±0.0005 


-0.7493 ± 0.0005 


-0.300 


-0.7635 ±0.0010 


-0.7182 ±0.0010 


-0.6846 ±0.0010 


-0.200 


-0.6702 ±0.0010 


-0.6313 ±0.0010 


-0.6032 ±0.0010 


-0.150 


-0.6112 ±0.0015 


-0.5764 ±0.0010 


-0.5513 ±0.0010 


-0.125 


-0.5764 ±0.0015 


-0.5445 ±0.0015 


-0.5213 ±0.0010 


-0.100 


-0.5366 ±0.0020 


-0.5077 ±0.0020 


-0.4865 ± 0.0024 


-0.075 


-0.4902 ± 0.0024 


-0.4641 ± 0.0020 


-0.4459 ± 0.0024 


-0.050 


-0.4299 ± 0.0026 


-0.4092 ± 0.0030 


-0.3934 ±0.0024 


0.050 


0.4273 ±0.0026 


0.4062 ±0.0030 


0.3915 ±0.0024 


0.075 


0.4877 ±0.0024 


0.4621 ±0.0020 


0.4435 ± 0.0024 


0.100 


0.5348 ±0.0020 


0.5055 ±0.0020 


0.4842 ± 0.0024 


0.125 


0.5749 ±0.0015 


0.5428 ±0.0015 


0.5201 ±0.0010 


0.150 


0.6093 ±0.0015 


0.5754 ±0.0010 


0.5505 ±0.0010 


0.200 


0.6693 ±0.0010 


0.6306 ±0.0010 


0.6023 ±0.0010 


0.300 


0.7625 ±0.0010 


0.7172 ±0.0010 


0.6838 ±0.0010 


0.400 


0.8367 ±0.0005 


0.7863 ±0.0005 


0.7488 ± 0.0005 


0.500 


0.8993 ±0.0005 


0.8445 ± 0.0003 


0.8035 ±0.0003 


0.600 


0.9540 ±0.0003 


0.8951 ±0.0003 


0.8513 ±0.0003 


0.700 


1.0024 ±0.0003 


0.9402 ±0.0003 


0.8937 ±0.0003 




= 33.3 ± 0.6 


Z^ = 26.0 ± 0.4 


Z^ = 21.3 ±0.3 




v B = (5.9±2.0)10- 5 


v B = (6.9 ± 1.4)10~ 4 


v B = (3.1±0.4)10" 3 




z l-l00 P = 32 g 


z i-ioo P = 26 3 
31 


z y oop = 21 g 



TABLES 

TABLE I. We report the values of 4>b(J) for the massless case as discussed in text. At the 
various (A , r Q ) we also show the results of the 2-parameter fit with Eq.(21) and the one loop 
prediction (25). 
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TABLE II. By using Eq.(21), we show the results of the 1-parameter fits for vb at A =0.8 and 
1.2 when is constrained to its one-loop value in Eq.(25). We also show the predictions from 
Eq.(24), (vb) T]cl , for yi = 2.07 ± 0.01 as determined from the fit to the data at A Q =1. 



A = 0.8 


A = 1.2 


r Q = -0.36 


r Q = -0.54 


= 32.90 = fixed 


Z<f, = 21.93 = fixed 


v B = (7.30±0.03)10" 5 


v B = (2.28±0.01)10~ 3 


x 2 _ 5 


X 2 _ 13 


d.o.f 21 


d.o.f 21 



(v B ) Th = (7.25 ± 0.03)10- 5 (v B ) Th = (2.29 ± 0.01)1Q- 3 
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TABLE III. We report the values of 4>b{J) as obtained by L. Cosmai and P. Cea with their 
16 4 lattice at A G = 1 and r Q = —0.45. Errors are statistical only. We also report the results of the 
fit with Eq.(29) and the one loop prediction (25). 



J 


A G = 1.0 r = -0.45 




0.100 


0.506086 ±0.91£- 


04 


0.125 


0.543089 ± 0.82E - 


04 


0.150 


0.575594 ± 0.82E - 


04 


0.200 


0.630715 ±0.62£ - 


04 


0.300 


0.717585 ±0.52£ - 


04 


0.400 


0.786503 ±0.UE- 


04 


0.500 


0.844473 ± 0A1E - 


04 


0.600 


0.894993 ± 0.39E - 


04 


0.700 


0.940074 ± 0.37E - 


04 


X 2 _ 3.2 
d.o.f 9-2 


= 26.323 ± 0.042 




v B = (5.783 ±0.013)10 


-4 



z l-loop = 26 _ 3lg 
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TABLE IV. We report the one-loop theoretical prediction J 1 1oop (</>b) based on Eq.(21), 
with the corresponding error AJ, which best fits the Monte Carlo data reported in Table III. 

a o ( rl-loop tEXP-,2 

= 26.319 = fixed and v B = 5.7943 10" 4 . x = ~ (aj? ■ 



4>B = 


0.940074 ± 0.37E 


-04 


J l - 


-loop 


= 0.699975 ± 0.82E 


-04 


jEXP 


= 0.700 


x 2 


= 0.09 


4>B = 


0.894993 ± 0.39£ 


-04 


J l - 


-loop 


= 0.600001 ±0.78£ 


-04 


jEXP 


= 0.600 


x 2 


= 0.00 


4>B = 


0.844473 ± 0A1E 


-04 


J 1 - 


-loop 


= 0.500042 ± 0.72E 


-04 


JEXP 


= 0.500 


x 2 


= 0.29 


4>B = 


0.786503 ± 0.44E 


-04 


J l - 


-loop 


= 0.400027 ±0.67£ 


-04 


jEXP 


= 0.400 


x 2 


= 0.17 


4>B = 


0.717585 ± 0.52£ 


-04 


J l - 


-loop 


= 0.299952 ± 0.65£ 


-04 


jEXP 


= 0.300 


x 2 


= 0.54 


4>B = 


0.630715 ± 0.62E 


-04 


J 1 


-loop 


= 0.199982 ±0.59^ 


-04 


jEXP 


= 0.200 


x 2 


= 0.09 


4>B = 


0.575594 ± 0.82E 


-04 


J 1 - 


-loop 


= 0.150010 ±0.64^ 


-04 


JEXP 


= 0.150 


x 2 


= 0.02 


4>B = 


0.543089 ± 0.82E 


-04 


J 1 - 


-loop 


= 0.124943 ±0.57£ 


-04 


jEXP 


= 0.125 


x 2 


= 0.99 



(f) B = 0.506086 ± 0ME - 04 ji-ioo P _ q. 100062 ± 0.54E — 04 J EXP = 0.100 \ 2 = 1-34 
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TABLE V. The same as in table IV for the leading-log fit. 



4>B = 


0.940074 ± 0.37E 


-04 


jLL 


= 0.696382 ± 0.82E - 


04 


jEXP 


= 0.700 


x 2 = 


1917 


4>B = 


0.894993 ± 0.39E 


-04 


jLL 


= 0.598116 ±0.78£- 


04 


jEXP 


= 0.600 


x 2 = 


580 


4>B = 


0.844473 ± 0.41E 


-04 


jLL 


= 0.499679 ± 0.72E - 


04 


jEXP 


= 0.500 


x 2 - 


= 19 


4>B = 


0.786503 ± O.UE 


-04 


JLL 


= 0.400981 ± 0.Q7E - 


04 


jEXP 


= 0.400 


x 2 = 


212 


4>B = 


0.717585 ± 0.52E 


-04 


jLL 


= 0.301938 ± 0.65E - 


04 


jEXP 


= 0.300 


x 2 = 


851 


4>B = 


0.630715 ± 0.62£ 


-04 


jLL 


= 0.202585 ± 0.60E - 


04 


jEXP 


= 0.200 


x 2 = 


1833 


4>B = 


0.575594 ± 0.82E 


-04 


jLL 


= 0.152692 ±0.65^ - 


04 


jEXP 


= 0.150 


x 2 = 


1675 


4>B = 


0.543089 ± 0.82E 


-04 


jLL 


= 0.127580 ±0.57£- 


04 


jEXP 


= 0.125 


x 2 = 


1910 



(f) B = 0.506086 ± OME - 04 J LL = 0.102582 ± 0.55.B - 04 J EXP = 0.100 \ 2 = 2175 
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TABLE VI. The values of pb{J) for the massless case are reported as discussed in the text. 
At the various values of A Q and r Q we also show the results of the 2-parameter fits with Eq.(29) 
and the one loop prediction (25). 



J 


A D = 1.0 r Q = -0.6 


A D = 1.5 r = -0.9 


A G = 2.0 r = -1.2 


0.050 


0.4110 ±0.0130 


0.3850 ±0.0116 


0.3753 ±0.0091 


0.075 


0.4686 ±0.0086 


0.4337 ±0.0086 


0.4176 ±0.0069 


0.100 


0.5133 ±0.0084 


0.4724 ± 0.0065 


0.4517 ±0.0055 


0.125 


0.5495 ± 0.0063 


0.5048 ± 0.0065 


0.4814 ±0.0054 


0.150 


0.5819 ±0.0060 


0.5332 ± 0.0053 


0.5063 ±0.0042 


0.200 


0.6377 ±0.0047 


0.5812 ±0.0040 


0.5500 ±0.0037 


0.250 


0.6844 ± 0.0039 


0.6217 ±0.0034 


0.5875 ±0.0029 


0.300 


0.7256 ±0.0030 


0.6571 ± 0.0026 


0.6190 ±0.0029 


0.350 


0.7612 ±0.0029 


0.6892 ± 0.0026 


0.6473 ±0.0026 


0.400 


0.7938 ±0.0025 


0.7178 ±0.0025 


0.6731 ±0.0025 


0.450 


0.8246 ± 0.0024 


0.7435 ± 0.0025 


0.6969 ±0.0023 


0.500 


0.8520 ±0.0023 


0.7683 ±0.0021 


0.7193 ±0.0023 


0.550 


0.8785 ±0.0021 


0.7911 ±0.0021 


0.7399 ±0.0018 


0.600 


0.9029 ±0.0019 


0.8124 ±0.0020 


0.7593 ±0.0017 


0.650 


0.9261 ±0.0019 


0.8330 ±0.0019 


0.7778 ±0.0017 


0.700 


0.9481 ±0.0018 


0.8520 ±0.0019 


0.7953 ±0.0017 


0.750 


0.9693 ±0.0018 


0.8703 ±0.0017 


0.8119 ±0.0015 


0.800 


0.9893 ±0.0016 


0.8879 ±0.0015 


0.8274 ±0.0013 




= 25.0 ± 1.6 


Z^ = 16.4 ±0.7 


Z^ = 12.3 ± 0.4 




v B = (1.4lJ;i)10- 3 


v B = (1.8±0.6)10~ 2 


v B = (5.6±0.8)10~ 2 




z l-loop = 2g 3 


z l-loop = n5 
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z l-loop = 13 A 



TABLE VII. By using Eq.(29), we show the results of the 1-parameter fits for vb at A =1.0 
and 2.0 when is constrained to its one-loop value in Eq.(25). We also show the predictions from 
Eq.(24), (vB) Th , for yi = 2.44 ± 0.03 as determined from the fit to the data at A G =1.5 . 



A D = 1.0 r Q = -0.6 


A 


= 2.0 r Q = -1.2 


= 26.32 = fixed 




= 13.16 = fixed 


v B = (7.09±0.11)10- 4 


v B = 


■■ (3.79±0.03)10- 2 


X 2 - 0.8 




X 2 - 7.3 


d.o.f 17 




d.o.f 17 



( VB )™ = (6.85 ± 0.09)10- 4 (^) Th = (3.89 ± 0.05)1Q- 2 
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TABLE VIII. We report the values of Pb{J) obtained on our 16 4 lattice for A D = 1.0 and 
r Q = —0.6. We also show the results of the 2-parameter fits with Eq.(29) and the one loop 
prediction (25). 

J X = 1.0 r D = -0.6 

0.050 0.4101 ±0.0098 

0.075 0.4659 ±0.0065 

0.100 0.5105 ±0.0052 

0.125 0.5470 ± 0.0046 

0.150 0.5774 ± 0.0035 

0.250 0.6808 ±0.0023 

0.300 0.7235 ±0.0022 

0.350 0.7601 ±0.0022 

0.400 0.7933 ±0.0022 

0.450 0.8237 ±0.0020 

0.500 0.8517 ±0.0019 

0.550 0.8779 ±0.0018 

0.600 0.9025 ±0.0018 

0.650 0.9256 ±0.0014 

0.700 0.9475 ±0.0013 

0.750 0.9681 ±0.0013 

0.800 0.9874 ±0.0012 
= 26.35 ± 1.33 
v B = (6.8i|;I)10- 4 

z l-loop = 
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